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ON SPIRAL MINIMAL SURFACES 



A. V. KISELEV AND V. I. VARLAMOV 

Abstract. A class of spiral minimal surfaces in E'^ is constructed using a symmetry 
reduction. The new surfaces are invariant with respect to the composition of rota- 
tion and dilatation. The solutions are obtained in closed form and their asymptotic 
behaviour is described. 



Introduction. In this paper we construct tlie two-dimensional minimal surfaces E C 
which are invariant with respect to the composition of the dilatation centered at the 
origin and a rotation of space around the origin. We present the solutions in closed form 
through the Legendre transformation, and we describe their asymptotic behaviour using 
a non-parametric representation. Also, we construct a special solution of the problem. 

The paper is organized as follows. In Sec. ^ we formulate the problem of a symme- 
try reduction for the minimal surface equation, and we express its general solution in 
parametric form using the Legendre transformation. In Sec. |2]we convert the reduction 
problem to the auxiliary Riccati equation and classify its solutions. Then in Sec. El we 
describe the phase portrait of a cubic-nonlinear ODE pP whose solutions determine the 
profiles of the minimal surfaces on some cylinder in E^. Finally, in Sec. 0]we construct 
two classes of the spiral minimal surfaces and indicate their asymptotic approximations. 
The shape of solutions that belong to the first class resembles the spiral galaxies. The 
second class of solutions is composed by helicoidal spiral surfaces with exponentially 
growing helice steps. 

1. The symmetry reduction problem 

Let ^, rj, and z be the Cartesian coordinates in space E'^. Consider the surfaces which 
are locally defined by graphs of functions, S = {z = 77) | (^, 77) G V C M^}. A 
two-dimensional surface E is minimal (that is, the minimum of area is realized by E) 
if the function x satisfies the minimal surface equation 

(1 + xj) Xaa - "^X^XvX^v + (1 + xl) Xvv = 0- (1) 

In the sequel, we investigate the properties of Eq. (|T)) and structures related with it up 
to its discrete symmetry C, ^ t], z —z. 

The Lie algebra of classical point symmetries for Eq. is generated PP by three 
translations along the respective coordinate axes, three rotations ^ j of the coordinate 
planes around the origin, and the dilatation A centered at the origin. Let \E'i be the 
(infinitesimal) rotation of the plane O^r/, and consider the generator = \E'i + A of a 
symmetry of Eq. (^Q). In jj, the problem of symmetry reduction for Eq. (Q) by the 
composition of rotation and dilatation was posed, although no attempts to find at 
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least one solution were performed. By constructions, each solution of the symmetry 
reduction problem is a minimal surface invariant with respect to the rotation of the 
plane O^t] and, simultaneously, the dilatation of the entire space E^. 

1.1. Solutions via the Legendre transform. The general solution of the reduction 
problem for Eq. by its symmetry is obtained |1] in parametric form using the 
Legendre transformation. Now we describe a two-parametric family of the 0-invariant 
minimal surfaces. Although, we claim that the reduction provides a special solution 
that is not incorporated in this family. The present paper is essentially devoted to the 
description of asymptotic properties of the two classes of these surfaces. 
First we recall that the generating section of the symmetry at hand is 

= X - - v)Xi - + v)Xv^ (2) 

see |21in] for details. The 0-invariance condition for xi^^v) is = 0. Next, consider the 
Legendre transformation 

C = {w = ^Xti + VXv -X, P = Xi, (1 = Xv}- (3) 

Now we act by the Legendre transformation C onto the system composed by Eq. ((T]) 
and the equation = 0. Then from we obtain the linear elliptic equation 

(1 + p^) Wpp + 2pq Wpg + (1 + q^) w^q = 0, (4a) 
and from the invariance condition = we get the equation 

w + qWp — pWg = 0. (4b) 

Let {g, ^) be the polar coordinates on the plane Opq such that p = g cos-i? and q = 
g sin {}. Then Eq. ()4b|l acquires the form w — -^w = 0, whence it follows that w = 
uj{g) ■ exp('i9). Substituting w{g, {}) in (jlaj), we obtain the equation 

g'-{l + g^) u;"{g) + gu;'{g) + u{g) = 0. (5) 

Its complex-valued solutions are 

^± = (2 + ^^)^ exp(^±arctan(-(l + ^^))"^±arctanh(-(l + ^2))^y (6) 

Their real and imaginary parts define a basis of real solutions for Eq. (jSj). 

The inverse Legendre transform is = = Wp, rj = Wg, x = Wp + IWq — w}. 
Rewriting in the polar coordinates g and we finally get 

duj 



n ^ • n' 

COS -i/H — sm V 



og 



■ exp {}. 



^ , ™. „ -exp??, 7]= — sin?? cost? -expt?, x 

og f? -I '-og ^ 

(7) 

Formulas ((Zj) provide a parametric representation of the generic 0-invariant minimal 
surfaces. Special minimal surfaces are defined by special solutions of Eq. (0) different 
from family (jUj). 

Remark 1. Representation (|7j) of the minimal surfaces based on the linearizing Legendre 
transformation, see Eq. Q, is not a unique way to describe open arcwise connected 
minimal surfaces E = = ^{p,q)} G E'^ in parametric form (e.g., we have p = X£. 
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and g = Xt? in formulas (jH)). An alternative is given through the Ennepert-Weierstrass 
representation |3] 

^ = {m I e = fo + Re $(A) d\, CeZcC, e E^}, (8a) 

here $(C) = {$i(C), '^'2(C)) "^3(0} is a triple of complex analytic functions that satisfy 
the constraint 

+ <l>^ + $2 = 0. (8b) 

We see that Re^ and Im(^ are the new parameters on the minimal surface S. It would 
be of interest to obtain the Weierstrass representation (jSJ for the 0-invariant minimal 
surfaces ([7j). 

1.2. Reduction of Eq. to an ODE. Let us return to the system composed by 
Eq. (P) and the constraint = 0. This symmetry reduction leads to one scalar second 
order ODE with a cubic nonlinearity. Namely, we have 

Proposition 1 (0). Let {p, 9) he the polar coordinates on the plane O^r] such that 
{z, p, 9) are the cylindric coordinates in space. Consider a minimal surface E G 
invariant w.r.t. symmetry Q of Eq. (^. Then E is defined by the formula 

z = p ■ h{9 — Inp), 

where the function h{q) of q = 9 — Inp satisfies the equation 

h" ■{h'' + 2) + h- 2h' - {h'f -{h- h'f = 0. (9) 

Corollary 2. Assume that E is a minimal surface invariant w.r.t. composition of 
rotation and dilatation, see (j21)- Let 11 = = h{q)\^_^^ he the profile defined hy E on 
the cylinder Q = {p = 1} C , here h is a solution of Eq. Q. Then the surface E 
is extended from Q along the intersections of the logarithmic spirals p = const ■ exp(^), 
^ e M, and the cones z = p ■ h{q). The (j)-invariant surface E has a singular point at 
the origin. 

The subsitution h = x, h' = y{x) maps Eq. (0) to the cubic- nonlinear equation 

(x=^ + 2)y^ = 2y-x + y^-ix- y)'; (10) 
equation (fTUjl defines the phase portrait of Eq. 0. 
Remark 2. Consider the family of equations 

(x' + 2)yy' = -x + ay + y^ + iy- xf (11) 

that contains (fTIH) as a particular case when a = 2. The change of variables = (x^ + 
2)y{x) in (fTTjl leads to the Abel equation u' = /3(x)m^ + /2(x)u^ + /i(x)u + /o(x), which 
is integrable in quadratures [Zj if a = 3. For a = 2 we have /o = — (^2^2)^ ' /i ~ 1^^' 
/2 = — (3x^ + 2), and fs = x^x"^ + l)(a;^ + 2). The Abel equation can be expressed in 
the normal form rj' = rf + H{x{^)), where r] = //(O? the prime denotes the derivative 
w.r.t. ^, and the function H{x) is constructed as follows. Let us make another change of 
coordinates, u{x) = uj{x)ri{^)-^, where uj{x) = exp(/(/i - S-)dx) = ^2/3(^2+1)5/6 and 



4 



A. V. KISELEV AND V. I. VARLAMOV 



^ = ^{x) = J ^i/3(-^2^^/)2/3 '^3^- The function x(^) is defined by inverting the expression 
for ^{x). Then the function H{x) is obtained from the relation 



whence we finally get 



/2^' /1/2 , 2/1 

I; 

3/2 



x^ + i y^^ -20_ 

_a;2 + 2y ■ 27x2(x2 + l)2(x2 + 2)2 

note that the inverse function x(^) must be substituted for x in ()12|1 . Hence we conclude 
that the coordinate transformation u{x) = ^2/3(^2+1)5/6 ■ viO + 3xix^+2}(l2+i) "^^P^ 
Abel equation to its normal form rj' = t]^ + H{x{C,)). 

1.3. Exact solutions of the limit analogue of Eq. ()10|) . In this subsection we 
consider the limit analogue of system (fTUI) whose right-hand side contains the terms 
that provide maximal contribution as i? = >y x"^ + y"^ 00 on Oxy. To this end, 
we note that for large R the right-hand side of Eq. (fTUj) is equivalent to its "cubic" 
component [y^ — (x — yY~\ ■ y~^ ■ x"^. Let us find exact solutions of the equation 

dy - {x- yf 



dx yx"^ 

Using the substitution s = Mn (fT!?jl . we obtain the equation 

ds 2s^ + 3s - 4s2 - 1 



(13) 



X 



dx s 

Note that 25^ - 45^ + 3s - 1 = (s - l)(2s2 - 2s + 1). Therefore we get 

- In — 5 = ln|(5x|, (14) 



(s-l)(2s2-2s + l) 2 2s2-2s + l 

whence we finally obtain the integral 2l^-2l+i ~ ^^^^ 6 >0. 

Suppose 5 = 0, then we have s = 1. Thus we get the solution y = x oi approxima- 
tion (fT^ . In Proposition ^1 below we prove that the diagonal is the asymptote for a 
solution of Eq. (fTUj). If 5 > 0, then we solve the quadratic equation w.r.t. s and obtain 

_ 1 - (5x2 ± 75 ■ |x| • ^1 _ ^3,2 
^ ~ 1 - 25x2 ■ 

Hence we conclude that 



m = ^^1^?^ ■ (15b) 

Vl-5a;2 + v^- |x| 

Solutions (fT3j) are defined for |x| < \x\ ^ Suppose x > 0, then we conclude 
that 

lim |/i(x) = +00. 

V2S 
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Finally, consider solution (jlSaj) with the Cauchy data {xo,yo), yo 7^ Xq, located on a 
large circle. Then the above reasonings yield that 



S 



I— — 1| 



We claim that the behaviour of solutions of initial equation (fTn|) is correlated with 
the solutions of limit equation (frS|) . see Remark |H1 on p. HU 

2. Solutions of the auxiliary Riccati equation 

In this section, we start to investigate the phase portrait of Eq. We show that 
the phase curves are described by solutions of the auxiliary Riccati equation (j221)- 
Using the topological method by Warzewski we conclude that Eq. has one 
unstable solution and a large class of stable solutions. On the phase plane for Eq. ^ 
these solutions are represented by a pair of centrally symmetric separatrices and by 
trajectories that repulse from the separatrices and have vertical asymptotes at infinity, 
respectively. This will be discussed in the next section. 

Consider the autonomous system associated with (^UJ, 

^ = 2y + yx^, ^ = 2y-x + y^-ix-yf. (16) 

Integral curves for Eq. (fTUI) are assigned to integral trajectories of system (fTHj) . Consider 
the linear change of variables x = 2{zi — Z2) , y = 2z\. Then system fllfj|l is transformed 
to 

^ = 21 +2:2 +4(2:5* + ^2^), '^ = -Zx^Z2^^Z2{2zl-ZxZ2^zl). (17) 

The linear approximation for (jl7p has an unstable focus at the origin; indeed, the 
eigenvalues are Ai^2 = 1 By a Lyapunov's theorem, nonlinear systems (jl6|) and (jl7|) 
will exhibit analogous behaviour near the origin. 

Next, we use the substitution z\ = rcosip, zi = r sin 99 in (fTTj). whence after some 
transformations we obtain the triangular system 

r = r + 4r^, (18a) 
= —1 + 4r^ sin (cosy? — sin (/?). (18b) 

For any Tq > the Cauchy problem for Eq. (jl8aj) has the solution 

/, X roexp(t) 
r(t,ro)= . (19) 



1 + 4r2(l -exp(2t)) 

whenever r(0, tq) = tq. Assume that r = ro>Oatt = for a solution {r{t), 
(p{t)}. From (fTU|) it follows that any solution of system (fTHj) achieves the infinity on the 
plane Oxy at the finite time t* = |ln(l + ^), and we see that t* 00 if ro +0. 
Obviously, we have — >■ — 1 for r — >■ +0, that is, the trajectories in a small neighbour- 
hood of the origin are the spirals that unroll clockwise. These reasonings describe the 
behaviour of the trajectories of systems (fTB|l and (fTTj). 
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Proposition 3. The trajectories Zi{t), Z2{t) are transversal to any circle centered at 
the origin. Therefore the phase curves {x{t), y(t)} for Eq. Q are transversal to any 
cetrally symmetric ellipsis that corresponds to a circle on the plane 0ziZ2 {the major axes 
of these ellipses are located along the diagonal y = x). Hence we deduce that system (tT6|l 
has no cycles and equation Q has no periodic solutions except zero, which corresponds 
to the stable point on the phase plane. A unique stable point (0, 0) of system (fTB|) is 
an unstable focus. The spiral-type phase curves for Eq. © unroll clockwise around this 
point. Any trajectory of system p6p that does not coincide with the origin achieves the 
infinity at a finite time. 

Remark 3. The equations in system p6|) do not change under the involution {x,y) 
(— —y). Hence for any trajectory of systems (fT^18|) its image under the central 
symmetry is also a trajectory. Further on, we investigate properties of the phase curves 
y{t)} up to this symmetry (or, equivalently, up to the transformation y9 i— > + tt 
of the angular coordinate ip). 

2.1. In this subsection we study the behavoiur of solutions of system ()18|) as r ^ +oo. 
Let us divide Eq. ()18b|l by Eq. p8a|) . Thus we obtain 

dip — 1 + 4r^ sin(y9(cos(/? — sinyj) 



dr r + 4r^ 

Substituting r = 4r^ in (pUj) . we get the equation 

dip — 1 + r sin(y9(cosv9 — sin(y9) 
~ 2r(l + r) ' 

After some trigonometric simplifications in the r.h.s. of Eq. (PT|) we arrive at 

d^ (V2-l)r-2-2x/2rsin2(y,-|) 



(20) 



(21) 



(EU) 



dr At{t + 1) 

Next, divide both sides in Eq. ()2m) by cos^(v9 — |) assuming that ip ^ + Tin, n G Z. 
Also, we put u = tan((y9 — |) by definition. Hence we finally obtain the Riccati equation 



^ = 6(r) - a(r)u2, (22a) 
dr 



where 



^ ' 4r(r + l) ' ^ ' 4r(r+l) ^ ' 

The Riccati equation (j^^ can be further transformed [TU] to the Schrodinger equation 
with the potential b/a and zero energy. 

2.2. The Warzewski theorem. First let us recall some definitions 

Definition 1. Let C i7 be an open subset of a domain Q. Consider the Cauchy 
problem 

y' = fit,y), (23a) 

y{to) = yo. (23b) 
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Suppose y = y(t) is a solution of (j^ . A point (to, l/o) ^ ^ H is called an exit point 
with respect to equation ()23a|) and the domain if for any solution y{t) satisfying ()23b|) 
there is a constant £ > such that (t, € for all t in the interval t^ — e < t < t^. 
An exit point (to^l/o) for the domain fi" is a sthct exit point if {t,y{t)) ^ whenever 
to < t < to + ^ for some £ > 0. We denote by ^2° the set of all exit points for the 
domain and let f2°g denote the set of all strict exit points. 

Now we analyze the behaviour of solutions of Eq. ()22j) as r ^ +oo. First we note that 
tan I = a/2 — 1; hence we put Ui = \/2 — l, U2 = l — \/2. The above notation corresponds 
to the angles Lfi = j and Lp2 = 0, respectively. Also, we set tq = 2(^/2 + 1) + So, 
where (5o > is arbitrary. 

Further, we introduce two closed domains Di and D2 in the right half-plane r > 
of the coordinate plane {t;u): we let 

Di = {{t;u)\t > To, 0<u<ui}, 

D2 = {{t;u)\t > To, U2<u< 0}. 



(24) 



Next, we note that the equality ^ = is valid on the curves u{t) = -L^^h{T) / a{T). 
Therefore we define the third domain 1^3 C -Di U D2 through 



^3= |(r,M)|r>2(V2 + l), -^h{T)/a{T)<u<^h{T)/a{T) 

It is easy to check that the inequality ^ > holds in the domain D3. 

By definition, put /(m, r) = 6(r) — a(r)M^. Let us describe the inclination field for 
Eq. on the lines u = ui, u = and u = U2. We have /(ui,r) = f{u2,T) = 

2r(r+i) ^ ^ whenever r > tq, and we also have /(O, r) = 6(r) > under the same 

assumption. This argument shows that all points on the lines u = and u = Ui are the 
strict entry points with respect to Di. Simultaneously, the lines -u = and u = U2 are 
composed by the strict exit points w.r.t. the closed domain D2. 

The Warzewski theorem and Example ^below are borrowed from Using them, we 
prove the existence of solutions of Eq. (j^^ that do not leave the domains Di and D2. 



Theorem 4 ([3 ). Let f{t,y) be a continuous function on an open set Q of points 
{t, y), and assume that solutions of system ()23ap are uniquely defined by the initial 
condition ()23b|) . Also, let G Q be an open subset such that = Further 
let S G be a nonempty subset such that 

• the intersection S nVl'^ is a retract of VP^, but 

• the intersection S* fl fi^ is not a retract of S . 

Then there is at least one point {to, yo) G 5* fl such that the graph of solution y{t) of 
the Cauchy problem is contained in on its maximal right interval of definition. 

Example 1 ([3]). Suppose that y is real and the function f{t,y) in system ()23a|) is 
continuous on the set VL which coincides with the whole plain (t, y). Let Vf* be the strip 
\y\ <h, — oo<t<oo. The boundary of is contained in Vt and consists of the two 
lines y = ±b. Assume f{t,b) > and /(t, -b) < such that = = 9^° n Q. 
Next, let 5* be the segment S = {it,y) | t = 0, \y\ < b}. Then S* fl fig consists of 
the two points (0, ±b); the intersection is a retract of of the set fig but is not a retract 
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of S. Theorem m yields the existence of a point {0,yo), \yo\ < b such that there is the 
solution of the Cauchy problem for system ()23a|) with the initial condition y{0) = y^. 
This solution satisfies the inequality \y{t)\ < b for all t > 0. 

From Example H that illustrates Theorem |3 we obtain 

Corollary 5. There is a solution u = iPi{t) of the Riccati equation (j22|) that does not 
leave the domain Di, which is defined in (f^ . for all r > Tq. Analogously, there is a 
solution u = ip2{'T) in D2 that does not leave D2. 

Remark 4. The Warzewski theorem guarantees the existence of a solution ipi in Di. We 
claim that there are infinitely many solutions of this class in our case. The asymptotic 
expansions of these solutions as r ^ 00 are specified in fl32|) . The convergence of the 
expansions can be rigorously proved, see Remark El on p. ^| Also, we claim that the 
solution iIj2 in the domain D2 is unique and unstable. 



Now we calculate the limits of solutions of the Riccati equation ()22p as r — >■ 00. 

Lemma 6. Let u{t) he a solution of Eq. (j22j) which is greater than u = V^i(t) for all 
T > Tq. Then we have 

lim(M(r) -^/'i(r)) = 0. (25) 

T— >00 

Proof. Assume the converse, inf (m(t) — iIji{t)) = Aq, where Aq > 0. Then the differ- 

T>TO 

ence w{t) = u{t) — '?/'i(r) of two solutions for Eq. satisfies the equation 

^ = -a(r)(^r) + ^i(r))-w;. (26) 
Integrating Eq. (|^. we obtain 

w(r,ro) = w(ro) ■exp(- / a{s){u{s) + i)i{s))ds. (27) 



TO 



Further recall that 

«(s)+^i(s) >u{s)-^i{s) > Ao 
whenever s >Tq and a(s) > 0. Therefore, 

-! a(s)(u(s)+Ms))ds < -Aq [ a(s)ds = -—((V2-l)\n^^ + 2\n 

J TO J TO 4 V ^0 + 1 



Consequently, we have lim w(r, Tq) = 0. This conclusion contradicts the assump- 

T— >+00 

tion Ao > 0. □ 

Similarly we prove that a solution u{t) tends to -u = ipiij) as r — > +00 if u{t) enters 
the domain Di and there it is located under the graph u = ipiir). 

Lemma 7. The following equality holds: 



T— >+oo r^+oo 



lim ipilr) = lim -i/6(r)/a(r) = -\/2 — 1 = Ui. 
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Proof. Suppose iPi{t) C DiCl for all r > Tq. Then the solution iPi{t) grows and is 
bounded. Therefore there is the limit 

lim ipii^r) = do. 

Next, let u = u{t) be a solution that enters the domain D^r\Di through the line u = Ui 
at a large ri and remains in it for all r > ri. If c/q < V2 — 1, then the limit of the 
solution u{t) is di > do, but the limits do and di must coincide by Consequently, 
do = di = V2-l. □ 

Lemma 8. The solution u = ip2{T) remains in the domain D2 \ for all r > Tq, and 
its limit at infinity is 

lim ^2(r) = 1-V2. 

Proof. If the solution ip2{^) enters the domain D2 fl -D3 at some ti > tq and remains 
there, then it grows and is bounded. Consequently, there is the limit lim ip2{^) = d,2, 

r— >+oo 

where d2 E {1 — v^; 0]. Hence for a large T2 and r > T2 we obtain 

^i[r) + V'2M > ^ • 

By definition, put w{t) = ipi^r) — ip2{T)- Then from Eq. it follows that 

whence we deduce 
w{t) = w{t2) ■ exp ^- J a{s){'ipi{s) + 'ijj2{s))ds^ < 

< w{T2)exp I - ^ ~ J ^ • / a{s)ds 1 . (28) 



V2 



Yet we see that the r.h.s. in ()28p tends to as r ^ +00. Therefore the l.h.s. in 
must also tend to zero which is impossible. Hence the graph V'2(t) is contained in 
D2 \ -D3 for large r, that is, the solution ?/'2 decreases and is bounded from below. This 
argument shows that 

lim V^2(r)= lim ( -./6(r)/a(r) ) = -{V2 - 1) = 1-^2. 

This completes the proof. □ 

Analogously, suppose the graph of a solution enters the domain and is located 
above the graph u = ip2{T)- Then it tends to the solution u = iPi{t). The proof is 
straightforward. A solution which is less than u = V'2(t) achieves —00 at a finite time. 
Using Eqs. and (j2Z|), it can be proved that the solution u = 'ip^iT) in the domain 
D2 \ -D3 is unique. The proof is by reduction ad absurdum. 
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2.3. Asymptotic expansions of the solutions iP{t). In this subsection we use the 
method of undetermined coefficients and obtain the asymptotic expansion in ^ for the 
solution (r) of Eq. that tends from above to Ui = tan | as r = ir'^ +00. Also, 
we get the expansion for ^2{t) that tends to U2 = — tan| at infinity. We emphasize 
that the expansion for ipl provides a solution different from iPi{t), which exists by 
Corollary El Indeed, the new function -ipl tends to the limit Ui = \/2 — 1 monotonously 
descreasing. 

Proposition 9. The solutions ipii'^) and ip2{'T) admit the following asymptotic expan- 
sions: 



^■;(r) = V2-l+ ''^'^-^> +0(l). (29a) 



feW-l-v^+ '^^'f -^> +0(l), (29b) 

The inequalities iI^\{t) > \pl — 1 and ip2iT) > 1 — \/2 hold for large r. 

Remark 5. Using the method of majorant series jHlE], we prove that bounded solutions 
of the Riccati equation are assigned to expansions (j^ . These solutions are real 
analytic in a neighbourhood of the infinity, and the radius of this neighbourhood can 
be estimated. 

2.4. The general case: u = ?/'(r). Now we construct the asymptotic expansions for 
all solutions of that tend to ±(-\/2 — 1). This time we use the expansions in r but 
not in r = 4r^. 

Let us re-write the Riccati equation (j^^ in the form 



du ^ 2(v/2-l)r2-l _ 2{V2 + iy + l , 

dr r^Ar^ + 1) r(4r2 + 1) ' " ' ^ ^ 

Its right-hand side is real analytic in Mf r > i. Suppose r > 1. Let us use the expansion 

Mr) =^0 + ^ + ^ + 0(1). (31) 
Substituting for u in (jHUIl . we get 

_w^_2w2_ _ V2-1 ,^ 1^1 x_ 

\_ /-, L_L^ \_{ 2 2WoWi w\^2WqW2 n 

2r ^ 4r2 ^ Vor^ • • -j + \ ^^2 ]^g^4 • ■ vy- 

Equating the coefficients of 1/r and l/r^, we obtain 

v^- 1 - (V2 + 1) -w^ = 0, wi = iV2 + 1) ■ wqWi. 

The former equation has two roots, wq,! = — 1 and Wo,2 = 1 — v^- First we let 
Wo = Wq^i] then the root of the second equation is an arbitrary real number! In this 
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case, we use the notation wi = C. Secondly, suppose Wq = 1^0,2; then a unique root of 
the second equation is Wi = 0. Now we equate the coefficients of ^, whence we get 

V2-I 1 V2 + I ^ , (V2 + IK 
-2W2 = - ^ [Wi + 2W0W2) + . 

If Wo = 72 - 1, then W2 = ^^^'^^ + ^ ■ C^. Alternatively, if Wq = 1 ~ V^, then 
W2 = y^^iy^zll^ This implies that for C = we obtain the two asymptotic expansions 
for the solutions ipKr) and ip2{T)i which depend on even powers of r and which were 
previously found in ()29p. Now we see that all other expansions involve r = ^ explicitly. 

By definition, put D = i [^/2{^/2 - 1) + (v^ + 1) ■ C^] . Then for any u{r) such that 
\mir^^oou{r) = ui we finally have 

uir,C) = V2-l + - + ^ + 0(\), CeR. (32) 

An analogue of Remark El is also true for ()32|) : using the method of majorant series 0, 
one readily proves the convergence of expansion ()32|) for large r to analytic solutions of 
equation (jSO)), and it is also possible to estimate the radius of convergence. 

Finally, we formulate the assertion about the behaviour of solutions of the Riccati 
equation (j221)- 

Theorem 10. //r G (0, 2[v^+l]), then all solutions of Eq. ()22|) decrease monotonously. 
Suppose T > tq. Then there is the unstable solution ip2{T) that tends from above to the 
limit U2 = — tan | as r — 00; the asymptotic expansion for ip2 'is given in ()29b|) . 

In the domain Di, see (j23); there are infinitely many growing solutions tpi{r) that 
tend to ui = tan | as r ^ 00. These solutions correspond to C < in expansion ()32|) . 

All solutions which are located between ip2{'r) arid tpi{r) tend to Ui as r = 4r^ 00. 
All solutions which are located under ip2 repulse from it; they decrease and achieve —00 
at a finite time. 

The solution ipUr), which is located above the domain Di, decreases and tends to ui 
as T ^ 00; its expansion is described by fl29aj) (or by formula ()32|) with C = 0). There 
are infinitely many other solutions '?/'i(r) which are greater than ipKr). These solutions 
also tend to Ui as t 00, and their expansions at infinity are given through (jH^ 
with C > 0. 

3. The phase plane Oxy 

In this section we describe the behaviour of phase curves for Eq. (jlUj) . Inverting the 
transformations introduced on p. Inland preserving the subscripts of respective functions, 
we pass from solutions iI>{t) of Eq. to solutions {x{t), yit)} of system (fTB|) . and next 
we obtain solutions y{x) of Eq. (fTUI) . For example, the solution ipKr) is transformed 
to yl{x), and the separatrix 1/2 is assigned to the unstable solution ip2- 

The exposition goes along the time t in ()16p : ffist we consider a neighbourhood of the 
origin, next we analyze the behaviour of the trajectories at a finite distance from the 
origin, and finally we describe their asymptotic expansions at the infinity (x^+y^ 00). 
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3.1. From Eq. ()18a|) it follows that system (fTBjl has an unstable focus at the origin. In 
a small neighbourhood of this point, all trajectories unroll clockwise. By Proposition 01 
each trajectory is transversal to the ellipses centered at the origin (we recall that these 
ellipses correspond to the circles r = const on the plane OZ1Z2)', the trajectories achieve 
the infinity at a finite time. 



Proposition 11. All extrema of the phase curves for system (fTI)|) are located on the 
straight line y = x/2. 

Proof. Solving the equation P-s^x, y) = {y — x)^ + y^ + 2y — x = with respect to y, we 
obtain the real root x/2. The quotient ^3(0;, y)/{y — x/2) = 2y'^ — 2xy + 2x^ + 2 has no 
real roots for y at any x. □ 



Consider the phase trajectories {x{t), vif)} of Eq. (fTUI) in a neighbourhood of the 
axis Ox. The parts of the trajectories near the points (x, 0) describe the graphs of 
solutions h{q) of Eq. (jU)) near their extrema. The trajectories are approximated with 
the circles of radius q{x) = \x\ ■ centered at the points (xq, 0), where xq = 
Obviously, we have Xq — > and ^(x) \x\ as |x| — > 00. 

Remark 6 (On inflection points). In the first and second quadrants of Oxy (and owing 
to the central symmetry of the phase portrait, in the third and fourth quadrants, 
respectively) there is the curve I that consists of the inflection points of the phase 
curves. The curve I is composed by two components Xi^2 that join in the first quadrant 
making a cusp; the point of their intersection is the nearest (w.r.t. the Euclidean metric 
on OZ1Z2) inflection point located on the spirals that unroll from the origin. 

The first component Ii consists of the inflection points where the trajectories start 
to repulse from the separatrix y2{x) (the unstable solution y2 tends from above to 
the ray y = x, see Proposition below) and turn towards their local maxima on 
the ray y = x/2. At the infinity x^ + — 00, the component Xi approaches the 
diagonal y = x from below. The second component X2 goes left from its intersection 
with Ti and contains the inflection point of the separatrix. Then the curve I2 enters 
the second quadrant. There, it describes the moments when the phase curves, having 
crossed the Ox axis at x < 0, repulse from the separatrix and turn upward, possessing 
the vertical asymptotes at infinity. 

3.2. Having achieved the inflection point, the trajectories of system approach 
the infinity following one of the three schemes below. The expansions for the solutions 
iPi{t), i^Kr), and ip2{T) of the Riccati equation (j221), which were obtained in the previous 
section, determine the asymptotes for the phase curves of all the three types. 

Proposition 12. The diagonal y = x is the asymptota of the separatrix y2{x) in a 
neighbourhood of infinity] the curve y = y2{x) approaches the diagonal from above. 
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Proof. Using expansion ()29b|) for the solution ^/'2(t) as r = 4r^ —>■ oo, and taking into 
account the transformation from Eq. (fTBj) to Eq. (|T7j). we obtain 

y — X = 2z2 = 2r sin ip = 2r sin((y9 ~ f ) + f ) = 

= 2rcosf cos((/?-f)- (tan(y?-f)+tanf) = 2r cos f (^2(r) + ^2 - l)/-\/l + ^i(r) = 

= 2-osf (^ll^ + 0(^))/Vl + ^i(r) ^ +0. 
This argument concludes the proof. □ 

The separatrix divides the trajectories in the first (consequently, in the third) quad- 
rant in two classes. The curves y\ of the first type repulse from the separatrix and go 
vertically upward, not leaving the quadrant. Other trajectories dive under the sepa- 
ratrix and achieve the extrema on the ray y = x/2. Then they intersect the axis Ox, 
whence they either turn down and have a vertical asymptote at — oo or reach the third 
quadrant. Then the whole situation is reproduced up to the central symmetry. 

Now we give a rigorous proof of these properties. 

Proposition 13. The coordinate axis x = is the vertical asymptote for the solution y\, 
which approaches it from the left. 

Proof. Using expansion ()29aj) and inverting the transformation Oxy ^ 0ziZ2, we let 
r — >■ CX3 and hence obtain 

X = 2{zi — Z2) = 2r{cosip — simp) = 

= -2V2rsin((/?- f) = -2v^r sin((^ - f ) - f ) = 
= -2V2cosf ■r(tan(¥?- |) - tan |) cos(y? - f ) = 

= -2v^cos f ■ r(V^i(r) - V2 + l)/^l+^j{T) = 

= -2V2cosl.r C^'^f^-'K o{^))/^l + rAr) ^ -0- (33) 

According to ()18aj) . the function r grows infinitely along the trajectories. Consequently, 
the function R{t) = \Jx'^{t) + y'^if) also tends to +00 along the solutions {x(t), y{t)}. 
By (jnni), the coordinate x tends to zero from the left, therefore the axis a; = is the 
vertical asymptote of the trajectory yl{x). □ 

Let us analyze the asymptotic behaviour of the trajectories yi{x), which correspond 
to the solutions tpi{r) of the Riccati equation (j22I)- Using expansion (jH^ for an arbi- 
trary solution ipi{r), we describe the asymptotes on the plane Oxy. We recall that the 
choice C > in formula (jH^ corresponds to solutions located above the graph ipKr). 
The expansion of the function ipl itself is given at C = 0, and the asymptote of its 
image on the plane Oxy was obtained in Proposition^] If C < 0, then the solutions 
■01 (r) are located between V'i(t) and 4'2{'t)- 

Consider the phase plane Oxy. The solutions of the class yi grow infinitely in the 
second quadrant if C > 0. Suppose C < 0, then the representatives of this class go to 
infinity between the axis Oy and the diagonal y = x in the first quadrant; recall that 
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the diagonal is the asymptote of the separatrix 2/2- Using an appropriate modification 
of the proof of Proposition [13 we obtain the estimate 

X = ^ (u{r, C) - tan = ^=^^= - + — + 0(— ) , 

where D = i [2 - ^2 + (v^ + 1) • C^] and expansion (jH^ is substituted for M(r, C). 
Hence we finally obtain 

Proposition 14. ^4// solutions of the class yi have vertical asymptotes, 

lim x{r, C) = -2V2C cos^ f = -( ^2 + 1)C. 

r— >oo 

r/ie graphs of solutions contained in the upper half -plane approach these asymptotes 
from the left. 

Remark 7. Without loss of generality, let C < 0. Then the vertical straight line 
X = — (-\/2 + 1) ■ C, which is located in the right half-plane, is the vertical asymptote for 
two trajectories of system (fTBjl . If a solution yi is greater than the separatrix 7/2, then 
it tends to the asymptote from the left. Another trajectory dives under the separatrix 
and approaches the same asymptote from the right in the fourth quadrant. The con- 
stants C > describe centrally symmetric curves in the third and second quadrants, 
respectively. 

Remark 8. From Propositions IT^HT^ it follows that for large R solutions of Eq. pO|) are 
approximated by the exact solutions of limit equation p3|) . The constants C in fl32p 
and S in dUD are related by the formula (a/2 + 1) ■ C = ±1 /y/26. 

4. The spiral minimal surfaces 

By construction, solutions h{q) of Eq. determine the profiles 11 = = h{q) \ 
p = 1, q E R} of the minimal surfaces E G on the cylinder Q = {p = 1}; we thus 
have n = E n Q. We further recall that a minimal surface in is extended from 
the profile 11 in agreement with Corollary |21 Yet, not each surface in assigned to a 
solution of Q is minimal. Let us study this aspect in more detail. 

4.1. The selection rule. In this subsection we formulate a rule that defines which 
components of the phase trajectories for Eq. (jU)) provide minimal surfaces. Using this 
rule, we conclude that the graphs {z = ^7)} of solutions x = P ' h{9 — Inp) attach 
nontrivially to each other. 

Let us recall a well-known property of the minimal surfaces [S]. A smooth two- 
dimensional surface in space is minimal iff at any point its mean curvature H vanishes. 
Hence each point of the surface is a saddle (of course, we assume that the surface at 
hand is not a plane). 

Further on, we denote by ^ the respective coordinate vector attached to a point of 
the cylinder Q. We see that the logarithmic spirals g = ^ — Inp = const are convex, and 
at any point of the surface the nonzero curvature vector has a positive projection onto 
^. Therefore we request that the curvature vectors at all points of the profiles H = 
{z = h{q), p = 1} C Q must have negative projections on ^. 
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The above condition upon solutions h{q) of Eq. © can be reformulated as follows: 
h' ■ h" > 0. This inequality is the rule for selecting the components of the phase curves 
on the plane Oxy. Namely, consider a neighbourhood Qr of a point (x, y) on a phase 
curve r. The minimal surface Er is assigned to this component of the curve if the 
coordinates x and y satisfy the inequality 

y^ ■ {y^ - (x - yf + 2y - x) > 0. 

Hence we obtain y > x/2 ioi y > and y < x/2 whenever y < 0, see the proof of 
Proposition [TTl If not all points of a phase curve satisfy this inequality (actually, each 
curve while it stays near the focus has to pass through the domains where it is not 
satisfied) , then several components of the profiles 11 and several minimal surfaces E are 
assigned to this curve. A continuous motion along the phase trajectory F corresponds 
to different (possibly, distant from each other) components of the profile 11. 

In Sec. El we described two classes of trajectories on the plane Oxy, the curves yi{x) 
with vertical asymptotes and the separatrix y2{x) with the diagonal asymptote y = x. 
Clearly, a small neighbourhood of the focus at the origin corresponds to oscillations 
of any profile, and some parts of the oscillating curves z = h{q) are prohibited by the 
selection rule. At infinity, the solutions of first type describe the finite size fragments of 
the profiles Hi = {z = h{q)}] these fragments have the vertical tangent at a finite value 
of the function h. Conversely, when the separatrix approaches the diagonal asymptote, 
it defines the exponentially growing profile 112 that turns infinitely many times around 
the cylinder Q. 

Example 2 (|1|). Let us choose the constant C < such that y*{ — {^/2 + l)C) = 0. 
Consider the following two components of the phase curves T = {y = yi{x, C)}: 

F_oo = {{x, y) I -{V2 + l)C > x{t) > 0, > y = yl{x) -oo}. 



{(x, y) I < x{t) < -(v^ + 1)C, y = yi(x) 



-oo}. 



-(y/2+l)C 

Also, let us use the focus (0,0), which corresponds to the trivial solution h = and 
hence amounts to the plane 2; = in E'^. Using these curves, we construct the closed 
profiles ni(C) C Q, see Fig. 1. 

h\ . h\ 



Fig.l 




Fig. 2 




By Corollary |21 the points R and S correspond to the logarighmic spirals on the 
plane z = in E^. Shifting the point R and the components ±cri towards S", we 
make the length \RS\ of the profile 11 comparable with 27r. Attaching the respective 
number of the profiles one after another on the cylinder Q, we place the edge R of 
each profile on the spiral S of the previous profile. Hence the resulting minimal surface 
becomes self-supporting. 

Example 3. Consider the following parts of the phase curves: a solution yi{x) as it 
tends to +00, another solution yi{x) of this class located helow the separatrix before 
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the extremum on y = x/2, and the separatrix y2{x) itself. Then we obtain the closed 
profile with finite cross-section in the upper half-plane h > 0. Indeed, we attach the 
corresponding curves ai, 5"i, and cr2 as shown on Fig. 2. 

4.2. The profiles U = {z = /ii,2(Q')} and their approximations. Now we obtain 
the asymptotic estimates for the profiles IIi. In this subsection we assume that the 
absolute value of the derivative h' = ^ (and, possibly, of the function h itself) is large. 
Let the constant C in (j32j) be arbitrary and put 

D = ^[2-V2 + C^-iV2 + l)] 

as before. Again, we use expansion (jH^ as r oo. Then we get the equivalence (~) 

h = X = 2{zi — Z2) = 2r(cos(/9 — simp) = 

2^/2 cos f ■r(V2- 1 -u(r,C)) . f ^ 



which holds up to O(^). The derivative h' is expressed through a solution of Eq. 
in the following way: 

h' = y = 2r cos ip = 2r cos | cos(v9 — |) (l — tan | tan(^ip — f )) = 

_ 2rcosf (1 - {V2-l)u{r,C)) 

As r 00, the above formula is equivalent (~) to 



2rcos=|(l-(v/2-l)'-(v/2-l)(^ + ^)) 



.(l + f)(2v/2-2)-(l + f)(v/2-l)(c+^l 



r / 



This reasoning shows that for large r we have 

Eliminating r from these formulas, we obtain the differential equation 

V2 h V2D 



h' 



2 c' + — 

I + V2 



By definition, put g = Then we get 
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Let US consider two cases: C = and C ^ 0. If C = 0, then we see that 

s/2{g^-{2-V2)) s/2 ^ n ■J , 

where A = const. Resolving this equahty w.r.t. g, we obtain 



±xl2-V2 + exp{-^^q + A). 

Now let C 7^ 0; then we express g{q) from Eq. (j^H). Therefore we geth= {V2 + l)g{q) 
in the form of an integral. 

Proposition 15. Suppose \h'\ — > (X). If C = 0, then the asymptotic behaviour of the 
solution h\{q) of Eq. Q, which is assigned to the phase curve y\{x), is described by the 
formula 



h\{q) ~ ±y 2 + \/2 + (v^ + 1)2 exp(V2(V2 - l)q + A), 

here A is an arbitrary constant. If C 7^ 0, then the approximation for h\{q) is given 
implicitly through the integral for Eq. (I34j) . 



\{y/2 - l)h + \C - Bf""^^ ■ \ {V2 - l)h + \C + Bf"" ^ = A-exp(2(2- y2)5g), 

here A>Q, B= sJ\C'^ + 2D, C e^, and D = \ [2 - ^2 + ■ {^2 + I)]. 

Finally, we use Proposition [T21 and describe the profile 112 assigned to the separa- 
trix ?/2- The asymptote y = x defines the equivalence h' ~ h as \h\, \h'\ 00. Conse- 
quently, the profile 112 is approximated by solutions of the differential equation h' = h. 

Proposition 16. Let h be large, then the function h2{q) grows exponentially m g G M: 
/i2 ~ A ■ exp(g), where A = const. 

4.3. Approximations of the minimal surfaces. From Propositions 1151 and fT6l we 

obtain the following assertion. 

Theorem 17. Let p, 6 be the polar coordinates on the plane O^rj and put q = 9 — In p, 

h = z\ . 
ip=i 

(1) Eor large values of h\q), there is a (p-invariant minimal surface El = {z = 
p ■ hl{6 — In p)} , which is approximated by the graph of the function 



/ (V2 + iyexp(V2(V2-l)e + A) , , , 

z = p.^2 + V2+'-^ —j^, ^ = --^- (35) 

Also, there is a class of the (p-invariant minimal surfaces Si that correspond to 
the asymptotic formulas for hi{q) in Proposition [T31 This class contains as 
a particular case. 

(2) Suppose that the absolute values I/12I are sufficiently large such that the phase 
curve h'2{h2) is near to the diagonal h' = h. Then there is the (p-invariant 
minimal surface S2 = {z = p- h2{6 — Inp)} whose approximation is given by the 
graph of the function 

z = Ap ■ exp(9 — Inp) = A ■ exp(^), A = const . (36) 
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The constant A G M determines the rotation of the surfaces around the axis Oz. The 
choice A < in (jHBj) defines the reflection symmetry z \—>- —z. Both surfaces E^, 
E2 C E'^ have a singular point at the origin. 

Finally, let us plot the graph of a spiral minimal surface determined by the profiles 11 
on Fig. 1, see Example El 



Fig. 3 

We conjecture that the spiral minimal surfaces may be relevant in astrophysics owing 
to their visual similarity with the spiral galaxy-like objects. Also, we indicate their ap- 
plication in fluid dynamics: these surfaces can realize the minimum of the free energy 
of a vortex. 
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